Qubit state monitoring by measurement of three complementary observables 
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We consider the evolution of a spin 1/2 (qubit) under the simultaneous continuous measurement of 
three non-commuting qubit operators a^, ay, a^. For identical ideal detectors the qubit state evolves 
by approaching a pure state with a random direction in the Bloch vector space and by undergoing 
locally isotropic diffusion in the perpendicular directions. The quantum state conditioned on the 
complete detector record is used to assess the fidelity of classically inspired estimates based on 
running time averages and discrete time bin detector outputs. 
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The needs of quantum computing/communication [1] 
are stimulating rapid progress in control of single quan- 
tum systems. Recent experiments demonstrate coherent 
manipulation of quantum systems, including Rabi oscil- 
lations and entangling operations with few qubits. An 
important direction for advanced quantum control is to 
realize continuous monitoring of a quantum system. The- 
ory of continuous quantum measurement [2t21 and ex- 
periments have been carried out on a number of 
systems. The quantum monitoring can be used topre- 
pare highly pure states and entangled states 0, 12-1^ 
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and for continuous error correction 

A particularly interesting case is when non-commuting 
variables are being measured simultaneously. In Ref. 
[l8| the signal cross-correlation for two such detectors 
of an evolving qubit was calculated. In Ref. [l^ Wei and 
Nazarov considered measurement outcomes for three de- 
tectors measuring a qubit in orthogonal directions. They 
analyzed the statistics of the integrated outcomes Vk for 
each detector and showed that if these outcomes happen 
to be sufficiently large, then the normalized vector v/\v\ 
is close to the Bloch vector of the actual qubit state. 

In this Letter we consider simultaneous continuous 
measurement of the qubit observables ax,(^y,^z, illus- 
trated in Fig. 1. The setup can be in principle realized 
with a trapped atom probed dispersively by optical cavity 



uz(t) 




FIG. 1: A qubit measured by three orthogonal detectors. 



fields. In contrast to [19|, we explicitly take into account 
the qubit evolution due to measurement, and analyze the 
problem of monitoring the qubit state using the measure- 
ment records. The three Pauli observables are comple- 
mentary, but the incremental changes of the quantum 
state due to the weak measurements carried out in in- 
finitesimal time intervals commute, and the simultaneous 
measurements contribute to purification of the quantum 
state three times faster than if only a single observable is 
measured. While the measurements drive the system to- 
wards the Bloch sphere surface (pure states), they cause 
locally isotropic diffusion in angular directions. 

If the observer has access only to detector read out 
signals integrated over finite time intervals, the back ac- 
tions associated with these accumulated signals do not 
commute, and the state of the qubit can only be approx- 
imately determined. The quality of the state estimate in 
this case is thus a measure of the role of complementarity 
of the observables detected. Comparing the exact qubit 
evolution with simple classically inspired ways of moni- 
toring, we show that running averages with an exponen- 
tial window can provide the fidelity of state monitoring 
up to 0.94. We also show that if the available measure- 
ment record is averaged over discrete time steps At, the 
monitoring fidelity decreases with At quite slowly. 

Model. We consider continuous measurement of the 
qubit observables ax, cfy, and by three linear detec- 
tors with output signals /^(t). Let Uk{t) = hit) — /o,fe, 
where /o,fe denotes the uniform average of the outcomes of 
detector k over the two qubit states, and let Auk denote 
the detector responses, i.e., the difference of the mean 
signal for the qubit states. We can then write 

Auk 



Uk{t) 



■Tr[p(t)(7fe] +^fe(t), k^x,y,z. (1) 
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The Tr[/5(7fe] expectation values, determined by the time 
dependent qubit density matrix p, are in the following 
represented as the Bloch vector r — {x,y,z). In Eq.((T]) 
^fe(i) are independent white noises with one-sided spec- 
tral densities Sk, {£,k^idt) — SkiSk/2. The qubit evolu- 
tion due to measurement by a linear detector (amplifier 
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with infinite gain) can be described by three parame- 
ters [1, [2^: the so-called measurement time ri„oas,fe = 
2S'fc/(AMfe)^, which determines the rate of quantum (in- 
formational) back-action, a factor K]^ describing the clas- 
sical back-action correlated with the output noise , and 
an ensemble dephasing rate Ffe, related to the single- 
qubit dephasing rate 7^ as Tu = jk + ^/2T^oiis,k+KlSk/4:. 
In this paper we are interested in the quantum back ac- 
tion due to measurements, and we assume the absence 
of classical back-action, = as well as the absence of 
any Hamitonian driving of the qubit. 

If the measurement is performed only by one a^- 
detector (F; = for I k), then the probability den- 
sity of its integrated result Ufc(r) = Uk(t)dt is 
Ptotiuk) = I]i Pii{^)Pi{uk), where the qubit density ma- 
trix p is written in the ak basis {i — 1,2) and Pi{uk) = 
yJrj'nSk exp{— [ufc -f (— l)*AMfc/2]^T/S'fe} are the Gaus- 
sian distributions for the basis states. Then the qubit 
evolution is given by the Bayesian quantum filter [6| 

Pu-M = A,(0)e-^''^(i-*-yP,(ufc)P,K)/Ptot(Sfc). 

(2) 

Qubit evolution with three detectors. In the case of 
three detectors measuring the qubit in the orthogonal 
bases, it is impossible to use the quantum Bayes rule for 
a finite r because the measurement back-actions do not 
commute with each other. Therefore we should apply 
Eq. ([2]) in the differential form (for small dt) in the three 
orthogonal bases corresponding to the measured observ- 
ables and then sum up the contributions to the qubit 
evolution. In this way we obtain the following equation 
in the Stratonovich form for the x-component of the qubit 
Bloch vector r{t) given the measurement record Uk{t): 

i = (1 - x'^){Aux/Sx)ux - xy{Auy/Sy)uy 

-xz{Auz/Sz)uz- {-^y +Jz)x. (3) 

Evolution equations for the components y and z can be 
obtained by cyclic permutation of variables in Eq. ([B]). 

Identical detectors. In what follows we consider the 
case of three identical detectors: Auk/Sk — Au/S — a 
(we assume a > 0) and 7^. =7 > 0. Then the qubit 
evolution Q can be rewritten in a vector form as 

r = —27 r + a {m (1 — r^) — [r X [r X It]] } , (4) 

where u = (Aii/2)r+^(t) is the vector of results, Eq. ([T|), 
and r = \r\. The evolution ([4]) is invariant under arbi- 
trary rotations and can be represented as evolution due 
to one-detector measurement along fluctuating random 
direction of it. It is interesting to note that while mea- 
surement of only single observable ak "attracts" the qubit 
state to one of the corresponding eigenvectors, the simul- 
taneous measurement of (t^;, CTj,, dz leads to no preferable 
direction in the Bloch space. 



The ensemble-averaged evolution is also isotropic: r = 
— 2rr, that is easier to see from the Ito form |21l] of (jj]): 

r^-2Tr + a{i{l-r'^)-[rx[rxi]]), (5) 

where F = 7 -I- Fq is the one-detector ensemble decoher- 
ence and Fq ~ {Au)^/AS — l/2rineas- We also introduce 
the efficiency (ideality) of the measurement rj = Fq/F. 

Transforming Eq. ([S]) to polar coordinates, we obtain 
the following evolution for the radial component r: 



f = 2Fo(l/r-r/r/)+a(l-r2)^„ 



(6) 



where (,r{t) ~ -^(i) is the noise component along r with 
the same spectral density: {(,r^rdt) — 5/2. In directions 
perpendicular to r Eq. ([5]) leads to a locally isotropic 
Brownian diffusion with coefficient a; correspondingly, 
the angular evolution of the projection onto the Bloch 
sphere surface has the diffusion coefficient a/r (this can 
be shown by using locally geodesic coordinates). 

In particular, for ideal measurement {r/ — 1) and 
pure initial state, the state remains pure [r = 1, 
see Eq. (|ni)] and the diffusion on the Bloch sphere 
can be described by the Fokker-Planck (FP) equation 
dp{9,(p)/dt = ToAg^tpp{9,ip) where Ag^^ is the angular 
part of the Laplacian. The solution of this equation 
at time r is p(e;T) = J2n=o ^e~"("+i)^/4p„(cos 9), 
where 9 is the angle from the initial state, Pn{z) are 
the Legendre polynomials, and V = 4Fot = 2r/Tincas is 
the variance. Obviously, for r 3> Tmoas the initial state 
is forgotten, and the distribution p(9;r) — )■ l/47r be- 
comes isotropic. We note that while the average state 
approaches the center of the Bloch sphere, the actual 
monitored qubit state remains pure, performing a ran- 
dom walk on the sphere. We also note that the back 
action quantified by the diffusion Fq decreases with in- 
creasing output noise S. 

Purification dynamics. As seen from Eq. ([5]), if r < 
77^/^, then on average f > 0, which means state purifi- 
cation. Uncertainty of our knowledge about the state is 
characterized by the linear entropy Sun — I — V, where 
V = 2Tr/5^ — 1 = is the state purity. For an ideal 
measurement, — I, and starting from a non-pure ini- 
tial state, the qubit will purify (T' — > 1) on a time scale of 
the order of Tmoas- For a non- ideal measurement, 77 < 1, 
purity will continue to fluctuate around a stationary av- 
erage value, {V)st < 1- 

To analyze the purification dynamics we use Eq. ^ 
to derive Ito equation for the purity: dV /dt = 2Fo[2(l — 
V/t^) + (1 - Vf] -f 2a(l - V)VV^r{t). The correspond- 
ing FP equation^ is ^ - ^ lA{V)p{r , t)] + 



\ -§p^ [B{V)p{r,t)] with coefficients A{r) = 2Fo[2(l - 
V/tj) + {1- V)^], B{V) = 8Fo7'(l - Vf, and initial dis- 
tribution p(7',0) — 6{V — Vo). At t ^ Tnicas the purity 
reaches a stationary distribution 



p,t{V,f])^N- 
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FIG. 2: Solid lines: averaged purity evolution {V)Fp{t) start- 
ing from the fully mixed state. Dotted lines: same evolution 
integrating Eq. (|8}. Dashed lines: evolution of (r) = {VP)fp- 
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FIG. 3: Monitoring fidelity ([9]) vs. window duration r for 
the rectangular (solid lines) and exponential windows (dashed 
lines) . Dotted lines show the uncorrelated product (r) (cos 0) t ■ 



where TV is the normalization. For 77 — > 1, pst{V,r]) ap- 
proaches the ^-function at = 1. 

In Fig. 2 we show the average purity (solid lines) 
{V)Fp{t) — Jq V p{V ^t) dV for measurement efficiencies 
77 — 1, 0.5, and 0.1, calculated by numerically solving 
the FP equation starting from the Bloch sphere center 
("Po = 0). The FP distribution p{V ,t) (not shown) has 
been also confirmed by the simulations using Eq. 

The purification dynamics can be approximated using 
the ensemble-averaged purification rate 13|, obtained 
from the above Ito equation, starting from a given purity: 



{dV)/dt = 2ro[2 (1 - Vh) + (1 - vf 



(8) 



For ideal detectors (77 = 1) this becomes 2ro(l — 'P){2> — 
V), that can also be easily obtained from the purifica- 
tion result [H m {dV)/dt = 2ro(l - V) {I - z^) for 
a z-detector by adding the contributions from measure- 
ments in X and y directions, so that 1 — — > 3 — r^. The 
purification by 3 detectors probing in mutually unbiased 
bases is on average 3 times faster than for a z-detector 
(averaging is over the directions in the Bloch space), since 
2-2 _ 7-2^3^ Also, the 3-detector purification rate {dV)/dt 
is isotropic, in contrast to the single-detector case, for 
which the adaptive measurement perpendicular to the 
spin vector leads to the fastest purification [l^, [l3| (Ref. 
[16| demonstrates that non-adaptive switching between 
random bases performs comparatively well) . It is impor- 
tant to note that the average purity {V)pp{t) differs from 
the naive integration of Eq. ^ (dotted lines in Fig. 2) 
because {dV) ^ d{V) and the purity distribution p{V,t) 
is generally different from (5- function. In particular, {V)st 
is slightly higher (for 77 7^ 0,1) than the stationary value 
(1 1/77) - y/{l + 1/77)2 - 3 derived from Eq. ©. 

Classically inspired state monitoring. Exact monitor- 
ing of the qubit state is realized by integrating the evolu- 
tion Eq. (j4]) given the measurement record Uk{t). How- 
ever, such real-time computation may be a challenge ex- 
perimentally, and therefore it is interesting to analyze the 
fidelity of simplified signal processing algorithms. To de- 



crease the noise component and reduce the bandwidth of 
signals given by Eq. ([1]), it is natural to average them over 
a running time-window: Uk{t) = g{t — t')uk{t') dt\ 
where g{At) is the window profile. We have considered 
(i) a rectangular window of duration r: (?(Ai) = for 
At < T and zero otherwise, and (ii) an exponential win- 
dow with decay time r: g{At) = exp(— Ai/r). The 
analyzed monitoring algorithm is very simple: at any t 
in the stationary regime (t ^ t, Tmeas) we estimate the 
qubit state as the pure state rcst(^) = The 
algorithm fidelity is defined as the time-averaged scalar 
product of this vector with the actual state r{t): 



F=2(Trpest/5)t-l- (r-ii/|M|)t 



(9) 



In Fig. 3 we show the fidelity F vs. the window du- 
ration r for the rectangular (solid lines) and exponen- 
tial (dashed lines) windows, calculated by simulating the 
evolution ([4]) for 77 = 1, 0.5, and 0.1. For 77 = 1 the 
fidelity reaches a maximum of -Fmax = 0.94 for the ex- 
ponential window with r — O.GTmoas (for the rectangu- 
lar window Fmax = 0.87 dX t — 0.9Tincas)- For small r 
the fidelity is suppressed due to large contribution from 
fiuctuations: \u\ ^ t~^/^, while for t > 77T,ncas it de- 
creases because signals from distant past loose their rel- 
evance to r{t). To analyze the latter effect quantita- 
tively, we have used Eq. ^ to find the signal-qubit cor- 
relations: {uz{t — At)z{t)) = (Au/2) exp(— At/77rinGas), 
{uz{t~ At)x{t)) = {uz{t — At)y{t)) ~ 0; other correlators 
are similar |24| . 

Since F = {r cos (f>)t, where 4> is the angle between u 
and r, the fidelity is bounded from above by the station- 
ary Bloch vector length (r) reached at t —J' 00 (Fig. 2, 
dashed lines). In Fig. 3 these bounds are shown as hori- 
zontal lines: (r) = 0.732 for 77 0.5 and (r) = 0.348 for 
77 = 0.1 (obviously, (r) = 1 for 77 = 1). It is interesting 
to see that with decreasing 77, the exponential-window 
Fmax approaches (7^). This means that at optimal r ei- 
ther u becomes practically aligned with r, (cos0)t — 1, 
or there is a significant correlation between fluctuations 
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FIG. 4: Monitoring fidelity vs. time step Ai for four algo- 
rithms processing time-discretized detector data (see text). 



of r(i) and (j){t). This correlation can be checked by com- 
paring F with the uncorrelated value (r)(cos0)t shown 
in Fig. 3 by dotted lines. We see practically no correla- 
tion at the optimal point for small rj, which means that 
(cos0)t — 1 (we have also checked this fact directly). 

Algorithms for time-discretized data. We consider now 
a situation when only the integrated detector signals 
= J*"_^^ui;{t') dt' are available at discrete time 
moments t„ = nAt, n = 0, 1, . . ., as in the experiment 
[lo| . Then exact monitoring of the qubit state is ob- 
viously impossible, especially if At ^ r^cas, and it is 
interesting to analyze performance of various monitor- 
ing algorithms in this case. Assuming ideal detectors to 
deal with pure states only, we have studied four algo- 
rithms, which output estimated states rest at moments 
t„, and calculated their fidelities (shown in Fig. 4) de- 
fined as F = (r • rest)n, similar to Eq. ([9]) but with av- 
eraging over moments tn- The actual evolution r(t) in 
this case is simulated via Eq. and we have checked 
that the fidelity does not depend on the inaccuracy of 
the initial state estimate. Algorithm 1 treats the vec- 
tor of measurement data it^"-* as a single measurement 
of the spin component along this vector and updates the 
qubit state using the quantum Bayes rule ([2]) in the cor- 
responding basis, changing at each time step. For small 
At the fidelity of this algorithm is F w 1 — 0.14Ai/T,„cas- 
Somewhat unexpectedly, even for At ~ T,„oas the fidelity 
is still quite good. Algorithm 2 at each step rotates the 
previous Bloch vector r(tn-i) towards the vector m^"^ by 
the angle A0 = (Au/S')u||"'' At, which is determined by 



the component m^"* of the vector tt'"' perpendicular to 
r(tn-i)- Even though for small At this is very similar to 
the Algorithm 1, the fidelity decreases more rapidly with 
increasing At. Algorithm 5 treats the three measurement 
outcomes as the results of sequential measurements of 
the three spin components and uses the Bayesian update 
rule accordingly. For Ai -> the fidelities of all three 
algorithms approach unity, though with different slopes. 
Algorithm 4 treats the data available at moments i„ in 



the same way as the running rectangular window (see 
the upper solid line in Fig. 3) and estimates the state 
as it("^/|tt(")|. Algorithm 4 suffers from large statistical 
errors for short At; however, for At « Tmeas the fideli- 
ties of algorithms 1 and 4 become practically equal, and 
for longer At the Algorithm 4 becomes the best among 
considered algorithms. 

In conclusion, state monitoring and purification by si- 
multaneous measurements of non-commuting observables 
has been described by quantum filtering theory. The 
shortcomings of simple, effective, algorithms reflect the 
difficulty of estimating quantum states from incomplete 
measurement data. The incompleteness of the time av- 
eraged or integrated data is due to complementarity and 
the non-commuting back action operations in the coarse 
grained limit of finite sampling times. Our analysis shows 
this very clearly and it quantifies the approach to perfect 
state estimation in the limit of continuous measurement 
and quantum filtering. 
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